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Abstract 

An analytic formula for the two-loop QCD correction of the decay H ^ 2'y \s 
presented. To evaluate all master integrals a 'Risch-like' algorithm was exploited. 
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1 Introduction 



The introduction of tiie Higgs particle into the Standard Model allows to describe the 
masses of all other particles by means of spontaneous symmetry breaking. The 'hidden' 
gauge symmetry also makes the model renormalizable. In this sense the Higgs particle 
is the crucial element in the SM of electroweak interactions. The mass of the Higgs is 
still unknown and thus the search for the Higgs is a tremendous experimental task. At 
present one deduces all possible properties of the Higgs from perturbation theory even 
in higher orders to be prepared for the interpretation of coming experiments at LHC. 
Among these are decay modes of the Higgs, which have been investigated in two loop 
approximation by many authors (for reviews see [1]). Our contribution to this problem 
is in the present paper a complete analytic calculation of the decay H ^ 2'y with the 
top quark as virtual loop. 

2 Results 

The partial decay width can be written in the form [1] , [2] 
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where r/ = Ml/Arnj, tw = Mjj/AMw- In the one loop order Aw{tw) is the dominating 
correction. We will focus our attention on At{Tt) and consider QCD corrections of order 
0{as) to this amplitude. It is convenient to write 
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The exact one-loop result is well-known [3] , [4] and reads 

^2 



Af' 



-A 



1 - 



4(1 - Ot)'- 



69, 



(3) 



where At = ^NcQf, v — 2 ^^^Gp^^^, Nq is the color factor, Qi is the electric charge 
of the t- quark and 
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The two-loop diagrams contributing to Afirt) are shown in Fig.l. There are six 
diagrams contributing to H^^ at the two-loop level (see Fig.l). 
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Fig.l. Two-loop graphs contributing to the QCD corrections to if 77 decay. 

The amphtude for the decay of the Higgs boson into two photons with polarization 
vectors e^(gi) and £,^(^2) has the following Lorentz structure: 



J2^t!i =Y1 (^t,i qiq2 9^'' + bt,i gi ^2 + ct,i gi ^2 ; 



(5) 



where Ct^i has no contribution for on-shell photons. In Eq. (5) the sum runs over all 
diagrams relevant for the decay H — > 77. Due to gauge invariance, we have J2i = 
— J2i bt,i- It is easy to find projectors for a^^j and bt f. 
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We calculated both, at^i and bt^i, in order to have an additional check for the correctness 
of our result. 

Tensor integrals and integrals with irreducible numerators were represented in terms 
of scalar integrals with shifted space-time dimension [5]. 

We introduced two auxiliary vectors ai and 02 allowing to obtain tensors with ki and 
k2 by differentiating with respect to ai,2, using the differential operator 
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where 



Q = — giai(ai(a3 04 + "5 + "e) + asae) — gia2(a3(tti + 02 + + aj) + 0:206) 
— g20.i(a2(a3 + 04 + ^5 + "e) + 04^6) — g20.2(a4(tti + 02 + «6 + "7) + «2a6) 
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Applying the integration by parts method, e.g., to integrals with selfenergy insertions 
in the fermion line 'scratches' corresponding fermion lines. For the remaining reduction 
of tensor integrals we apply (8) to the scalar integrals, setting Oj = after differentiation. 
It is also convenient to consider integrals of a more general type with the topology given 
in Fig. 2. To deal with six-line diagrams (Fig.l), the corresponding ctj are put to zero. 
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Fig. 2. Generic topology of the two-loop graphs 

With the help of recurrence relations one can reduce any integral with six lines to a 
combination of integrals with at least one line contracted. The two-loop correction in 
the on-shell renormalization scheme reads 



A 



(1) 



-At 



■et{i + et){i + ei) 



1 

48 



-A[uM) - Cs) InOt + 36Li4(-^t) + 27Li4( 



+ 



+ 2Li2(-0t)+C2j In' et- 
9t{Wl - 76^ + 2Wt + 3) 



+ 



+ 



2(1-^0' 



2eti39^ - 29t + 3) 
(1 - OtV 



12(1 - 9tY 
\netU2iet) 



9 



(1 - Otr 

9t{9^ - U9t + 1) 



C3 + 



391 



'l-9tY ' {l-9tY"^ - {\-9tf 

The expansion at = agrees with the known result [6], [7] 
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and the asymptotic behavior near threshold ~ 1 is: 
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At > 4mf the imaginary part of (10) is in agreement with the result of [8]. 
The only master integral of the transcendentality four contributing to the result is: 
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where P/T^ = — + ie and qf = q2 = 0, (gi — ^2)^ = Mjj. For the calculation of this 
result a method to some extent resembling the 'Risch algorithm' [9] was developed. To 
conclude, we stress that the availability of analytic results is of great value since they 
allow to cover the whole kinematical domain with one expression, serve as checks for 
approximation, allow to represent the analytic properties of the amplitudes and can also 
be used for numerical evaluation in all domains of interest. 
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